x=load(argv{1});
y=load(argv{2});
targets=x(:,1);
pos=sum((targets+1)/2);
neg=length(targets)-pos;
tp=(pos+1.0)/(pos+2.0);
tn=(1.0)/(neg+2.0);
for i=(1:length(targets))
    if (targets(i)<=0) targets(i)=tn;
    else targets(i)=tp;
    end
end
[t,b]=logistic_regression(targets,x(:,2));
z1=1./(1.+exp(t-b.*y(:,2)));
[t,b]=logistic_regression(targets,x(:,3));
z2=1./(1.+exp(t-b.*y(:,3)));
[(y(:,1)+1)/2,z1,z2]
